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ABSTRACT 

Aims. We investigate the physical properties and spatial distribution of Carbon Monoxide (CO) gas in the disks around the Herbig 
Ae/Be stars HD 97048 and HD 100546. 

Methods. Using high-spectral-resolution 4.588-4.715 pm spectra containing fundamental CO emission taken with CRIRES on the 
VLT, we probe the circumstellar gas and model the kinematics of the emission lines. By using spectro-astrometry on the spatially 
resolved targets, we constrain the physical size of the emitting regions in the disks. 

Results. We resolve, spectrally and spatially, the emission of the 13 CO v(l-0) vibrational band and the 12 CO v = 1 - 0, v = 2 — 1, v = 
3-2 and v = 4 - 3 vibrational bands in both targets, as well as the 12 CO v = 5 — 4 band in HD 100546. Modeling of the CO emission 
with a homogeneous disk in Keplerian motion, yields a best fit with an inner and outer radius of the CO emitting region of 1 1 and > 
100 AU for HD 97048. HD 100546 is not fit well with our model, but we derive a lower limit on the inner radius of 8 AU. The fact 
that gaseous [OI] emission was previously detected in both targets at significantly smaller radii suggests that CO may be effectively 
destroyed at small radii in the surface layers of these disks 

Key words, circumstellar matter, stars - pre-main-sequence, stars - HD 97048, stars - HD 100546, protoplanetary disks 



1. Introduction 

We present a study of warm carbon monoxide (CO) gas in 
two circumstellar disks surrounding youn g intermedia t e-mas s 
stars. CO is commonly detected in disks (Naii ta et all [2003), 
and often used as a tracer of the physical conditions in the 
disk and the kinematics of the gas. Ro-vibrational transitions of 
warm (T - 120-101 K) CO gas have already been found by 
iMitchell et al.l (fT990|) in 8 ou t of 9 surveyed embedded infrared 
sources. iBrittain et al.l (120071) had a 100 % detection rate for the 
ro-vibrational CO transitions in the disks around 9 intermedi- 
ate mass stars with optically thick inner disks, while their detec- 
tion rate is only 1 out of 5 for disks with an o ptically thin inner 
disk, a case they explain by UV fluorescence. IB lake & Boogerj 
(2004) ha ve detect e d CO gas in 5/5 Herbig Ae stars and recently, 
iPontoppidan et al.l d2008l) have spatially and spectrally resolved 
the 4.7 /im CO ro vibrational lines in the disks around 3 young 
solar mass stars with known dust gaps or inner holes, and de- 
tected CO well inside the dust gaps in all disks. 

HD 97048 (Sp. type AOpshe, (Ivan den Ancker et all [19981) ) 
is a well-studied nearby Herbig Ae star, whose spectral energy 
distribution (SED) displays a large infrared excess above the 



emission from the star, commonly explained by emission from 
warm dust in a circumstellar disk. The dis k is a flaring disk ac - 
cording to the clas sification proposed by Mee us etal] d200ll) . 
lLagage etai] (120061) resolved the disk around HD 97048 in the 
mid-IR and determined the flaring angle. PA H emission from 
radii u p to 200-300 AU has been detected by Ivan Boekel et aTJ 
(2004). Ha bart et al.l d2004l) have resolved emission from the 
3.43 /im and 3.54 //m diamond bands from within the inner 15 
AU in the disk. Circumstellar gas has also been detected, e.g. by 
lAcke & van den Ancketj d2006l) who detect spatially and spec- 
trally resolved [OI] emission from a few tenths to several tens 
of AU. HD 97048 is one of the only two He rbig stars where the 
H2 S( l) line at 17.035 /im has been detected (Ma rtin-Zaidi et all 
2007), coming from warm (T ex < 570 K) gas within 35 AU. 
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HD 100546 (Sp. type B9Vne dvan den Ancker et all [19981) ) 
also has a SED characteristic for a flared disk. Its IR spectrum in- 
dicates ajrigh_cry_stalh^e sihcate dust fraction in the disk surface 
layers dBouwman et all [20031) . There is observational evidence 
for a disk gap, or 'second wall', at » 10 AU from the centra l star. 
This was derived based on the SED (Bouwman et al., 2003), and 
confirmed by nulling interferometry dLiu et al.Ll2007l) . The [OI] 
6300 A line was studied by lAcke & van den Ancketi ((2006), 
who also find evidence for a disk gap at that distance, possi- 
bly induced by a very low mass stellar companion or planet of 
« 20 Mj up , or a planetary system. The latter authors spatially and 
spectrally resolve the [OI] emission, and determine it to orig inate 
from a few tenths to several tens of AU. iGeers et ail (|2007) have 
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Table 1. Astrophysical parameters of the programme stars 



HD 


log r efl 


log L bol 


M 


Dist. 




PA 




log [K] 


log [La] 


Wq] 


[pc] 


n 


[°] 


97048 


4.00 


1.42 


2.5 ±a2 


180^ 


42.8 + ""' 6 
42 ±» 


160± i9, c 


100546 


4.02 


1.62 


2A ±o.i 


103~ 


145 ±5 - d 



igooi) 



° Ivan den Ancker et al.l d 1 9981) lLagage et alj 

lAcke & van den AnckeJ d2006l) lArdila et al.l d2007» 

Table 2. Log of spectroscopic observations at June 16 th 2007 



Object 


UT start 


Ref. wl. 


N exp 


Exp. time 


PA range 




hh:mm:ss 


[nm] 




[s] 


[°] 


HD 97048 


00:41:15 


4662.1 


20 


15 


41.1-45.6 




00:58:43 


4676.1 


20 


15 


46.0 - 50.3 


HIP 052419 


02:36:14 


4662.1 


2 


15 


85.1 - 85.3 




02:38:47 


4676.1 


2 


15 


85.7 - 85.9 


HD 100546 


02:58:10 


4662.1 


6 


15 


74.8 - 76.3 




03:05:44 


4676.1 


6 


15 


76.7 - 78.2 


HIP 061585 


03:49:57 


4662.1 


2 


15 


72.5 - 72.7 




03:52:30 


4676.1 


2 


15 


73.1 -73.3 



detected PAH emission from 12 + 3 AU, and the circumstellar 
material has been spatially resolved many times, revealing rich 
structures in the disk. 

In this paper, we present high spectral resolution observa- 
tions with milli-arcsecond spatial precision of CO emission orig- 
inating from the disks around HD 97048 and HD 100546. 

2. Observations and data reduction 

High-spectral-resolution (R « 94000, determined from telluric 
lines) spectra of HD 97048 and HD 100546 were obtained on the 
early morning of June 16' /! 2007 with the VLT cryo genic high- 
resolu tion infrared echelle spectrograph (CRIRES0, Kaufl et al. 
(2004)). Adaptive Optics were used to optimize the signal-to- 
noise and the spatial resolution of the observations. These cover 
the wavelength range of 4.588 - 4.711 /mi, and contain many 
fundamental ro-vibrational transitions of the CO molecule. The 
observations were made with a slit width of 0.2", and with 
the slit rotated along the parallactic angle. HD 97048 and HD 
100546 were observed with 2 wavelength settings, each for 10 
and 3 minutes respectively, and the slit Position Angle (PA) dur- 
ing the observations spans a range of 41.1 - 50.3 degrees for HD 
97048 and 74.8 - 78.2 degrees for HD 100546. This spread in 
PA, close to the minor axis of the disks, is not optimal for our 
astrometry purpose, but was enforced by the observational cir- 
cumstances. The data were reduced using the CRIRES pipeline 
V1.5.C0, which performs wavelength calibration, background 
subtraction and Hatfield correction. 

To correct for telluric absorption we have observed a telluric 
standard star directly after each science target: HIP 052419 (Sp. 
type B0V) for HD 97048 and HIP 061585 (Sp. type B2IV) for 
HD 100546. Each standard was chosen to be as close as possi- 
ble on the sky so that their spectra are affected by similar atmo- 
spheric condit ions. We compared their spectra to an appropriate 
iKuruczld 19911) stellar atmosphere model to determine the instru- 
mental response. The optical depth of the telluric lines of the 
standard was scaled to that of the science target. The science 
spectra were derived by division with the instrumental response. 



1 http://www.eso.org/sci/facilities/paranal/instruments/crires/ 

2 http://www.eso.org/sci/data-processing/software/pipelines/index.html 



A sub-pixel wavelength shift was manually added to reduce tel- 
luric residuals (i.e. spikes due to minor wavelength mismatches). 
Some telluric absorption lines are fully saturated, causing prob- 
lems with the division of the spectra (i.e. division by 0). In the 
further data reduction, we will disregard these areas. We have 
shifted the spectrum to center the v = 2 — 1 lines on their rest 
wavelength. The spectra are shown in Fig.Q] 



3. Observational Analysis 

To check whether the targets are spatially resolved, we per- 
form spectro-astrometry on the reduced 2D frames. We deter- 
mine the location of the photocenter of the PSF and the spatial 
FWHM as function of wavelength (or equivalently velocity for 
each line), by fitting a Gaussian to the spatial profile in each 
spectral column. The relative position of the photocenter to the 
continuum, the Spatial Peak Position (SPP), is a measure for 
the spatial offset of the emitting region compared to the con- 
tinuum. Our observations are made with the slit positioned over 
a sma ll range in position angles. We note that Pontop pidan et al.l 
(2008) have observed three nearby T Tauri stars at 6 different 
PA's with CRIRES and the same observational set-up as em- 
ployed by us, and find, after pairwise subtracting parallel and 
anti parallel spectro-astrometric signals, that no telluric feature 
produces false spectro-astrometric signals in excess of 0.2-0.5 
mas. Since we do not possess anti-parallel spectra, we show the 
telluric absorption in Figs.|2]and[3] and note, by comparing with 
CO lines un-affected by telluric absorption, that their effect on 
the astrometric signal is small. We compare the SPP in the CO 
lines to the SPP of the underlying continuum (SPP=0 by defi- 
nition). If the target is resolved, and the photocenter of the CO 
emission region is offset to that of the dust continuum emission 
from the disk, the SPP is expected to be non-zero. The photo- 
center displacement can be derived from the measured SPP if 
we correct for dilution by the continuum flux, by multiplying 
the SPP with a factor 1 + F c {v)/F,(v ) with F c (v) and Fj(v) th e 
continuum and line flux at velocity v dPontoppidan et all 12008). 

The Gaussian fit to the PSF also yields a spatial FWHM. If 
the continuum is unresolved, this quantity is a measure for the 
achieved effective spatial resolution in the direction parallel to 
the slit. When resolved, the FWHM allows us to determine the 
spatial extent of the emission. Outside the spectral lines it then is 
a direct probe for the size of the continuum emitting disk, whilst 
inside the lines it is a function of both the line- and continuum 
emission, and its value a lower limit to the CO emission. Thus, 
measuring the SPP and the FWHM of the PSF in the continuum 
and the CO lines allows us to determine [ 1 ] the size of the con- 
tinuum emitting region, and [2] the spatial offset of each CO line 
with respect to the continuum. This procedure is shown in Figs. 
[2]and[3] 



4. Results 

We detect in the spectrum of HD 97048 the 13 CO v = 1 - 0, as 
well as in the 12 CO v = l-0tov = 4- 3 vibrational bands. All 
lines are spectrally resolved with a remarkably constant FWHM 
of 15 km s . The Half Width at Zero Intensity (HWZI) is more 
difficult to estimate due to the continuum S/N, but we make a 
conservative estimate of 15 km s corresponding to an inner ra- 
dius of 4.5 AU, assuming the stellar parameters and inclination 
as given in Table [TJ Average line profiles and astrometric signals 
of the 12 CO v = l-0tov = 3- 2are shown in Fig. [5] and 
Table[3] The 4.6 /urn continuum emission of the disk around HD 
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Fig. 1. Reduced CRIRES spectra of HD 97048 (top) and HD 100546 (bottom). The black spectrum has been used for the line 
measurements, whereas the red line shows the areas that have been ignored due to large errors, mostly caused by the correction 
for saturated telluric lines. The line identifications plotted above show detections of the isotope ( 12 C0 and 13 CO), the vibrational 
transition (e.g. v = 1-0), and the rotational transition (R(xx) or P(xx)). The broad hydrogen recombination lines Pfund B at 4654 
nm and Humphreys e at 4673 nm are also shown. 



97048 is resolved. Its FWHM of 0.231" + 0.004" corresponds 
to 23 ± 3 AU, after correcting for the FWHM of the PSF. The 
latter (0.192" + 0.001") was estimated from 10 unresolved tel- 
luric standards, all sufficiently bright for the AO loop to close at 
maximum frequency, observed with the same instrumental set- 
tings, and within 2 days of our observations. The red- and blue- 
shifted wings of the SPP form a sinusoidal signal, typical for an 
inclined circumstellar disk. 

In the spectrum of HD 100546 we detect the 13 CO v = 1 - 
0, as well as in the 12 CO v = 1 - to v = 5 - 4 vibrational 
bands. All lines are spectrally resolved and have a FWHM of 
15 km s' 1 , irrespective of the transition. A conservative estimate 
for the HWZI is 17.5 km s _1 corresponding to an inner radius 
of 3.3 AU. The average line profiles and astrometric signals of 
the 12 CO v - l-0tov = 4 — 3 transitions are shown in Fig. [7] 
and their values in Table [3] The 4.6 jum continuum emission of 
the disk around HD 100546 is resolved. Its FWHM of 0.218" + 
0.002" corresponds to 1 1 + 1 AU, after correcting for the FWHM 
of the PSF. The SPP is very asymmetric, i.e the negative shift of 
the blue wing in the SPP is larger than the positive SPP of the 
red wing. This asymmetry is also seen in the FWHM and is more 
pronounced for the lower vibrational transitions. 

5. A Simple Model 

A powerful method to determine the extent of the emitting gas is 
to model the kinematics of the line- and astrometric profiles. We 
construct the average profiles for the v = 1— 0tov=4-3 vibra- 
tional bands in Figs.|5]and|7] and fit these to the model shown in 
Fig. |4] In this model, the gas is in Keplerian orbit in a flat disk 
with known inclination and PA around a star with known stel- 
lar mass and distance. Further, the intensity of the emission de- 
creases as I(R) — I m , with I[ n the intensity at the inner ra- 
dius R m , and R the radial distance from the star dCarmona et all 
|2007|) . The astrometric signal is consequently determined by fit- 
ting a Gaussian to the simulated PSF. Because the stellar pa- 
rameters are all well constrained, the free parameters are a, R m , 
and R out . a describes the change in intensity of the CO emission 
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Fig. 2. In all three windows the thick black line represents the 
spectrum of HD 97048 (orange thin line) smoothed over 2 pix- 
els, and the red line denotes the same data for the telluric stan- 
dard. Shown with the vertical dotted lines are the locations of 
the (from left to right) 12 CO v = 4 - 3fl(23), v = 3 - 2R(U), 
v = 2 - lfl(06), v = 5 - 4fl(33) and 13 CO v = 1 - 0R{\2) emis- 
sion lines. Top: Close up of the spectrum of HD 97048 (as shown 
in Fig. Q). Middle: The SPP of HD 97048 and of the telluric 
standard, shifted down by 10 mas. Bottom: The spatial FWHM 
of the PSF for each spectral row for HD 97048 and the telluric 
standard. 



as function of radius, where a — 2 reflects an intensity directly 
proportional to the stellar radiation field assuming the surface 
density has no radial dependence. We compare the line profile, 
the SPP and the FWHM with their simulated counterparts for the 
parameter-space 1AU < R m < 29 AU, 25 AU < R out < 200 AU 
and < a < 5.5, and quantify the best fit by minimizing the 
reduced^- 2 . We estimate the variance on the data from the con- 
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Fig. 3. Same as figure|2]for HD 100546 

Table 3. Spatial information for the averaged vibrational tran- 
sitions shown in Figs [5] and [7] and the radii used to model the 
kinematics. 



Fig. 4. Projected images of the simulated CO with the best fit 
parameters for HD 97048 and HD 100546, together with the slit 
of the spectrograph at the time of observation. North is up, East 
is left. 



The spatial distribution of gas in HD 100546 is more com- 
plex than the 'simple' homogeneous disk in Keplerian motion 
we model, and the best fit values for each transition vary, with 
a lowest value of > 4 AU and a > 1.5. If we assume that all 
vibrational states have the same spatial distribution and origin, 
and we only fit the line wings, a more likely value of R m = 8 
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" Note that this FWHM is a function of continuum + line emission, 
but that the SPP has been corrected for the continuum dilution, and thus 
is a direct function of the CO emission. All values are directly derived 
from the 2D spectrum, and not the maximum value resulting from an 
alignment of the slit with the semi-major axis. 



tinuum adjacent to the line, and determine the best fit from the 
v = 2 — 1 line profile, because the v = 1 — data is affected by 
telluric absorption, and because this transitions has a higher S/N 
compared to the higher vibrational transitions. Since the con- 
tinuum corrected SPP is similar for each vibrational transition 
within the separate stars, we assume that all vibrational transi- 
tions have the same spatial origin. For HD 100546, and in lesser 
extent HD 97048, our model produces a more distinctly double 
peaked line profile than the observed profiles. This phenomenon 
is also observed for 9 out of 60 protoplanetary disks (Bast et al., 
in prep.), and suggests there is an extra low projected velocity 
gas component present in those disks. We have repeated the fit- 
ting procedure for only the part of the line profile sensitive to the 
inner part of the disk, the line wings (| Vp ro j | > A.5km s -1 ), and 
find no difference for the best fit. 

The best fit to the v = 1-0, v = 2 — 1, and v = 3 - 2 vi- 
brational transitions of HD 97048, are R& = 11 AU, R out > 100 
AU and a = 3.00. We show the surfaces for the v = 2 — 1 
fit to the line profile, the SPP and the FWHM in Fig. |6] These 
values have been derived using a PA of 180°, estimat ed from di 



to the p resence of a disk gap a nd consequent "second wall" at 
TO AU dBouwman et al.L 120031) . The average profiles of all vi- 
brational bands have a FWHM of 15 km s' 1 , hinting that they 
share a common spatial origin. The [OI] emission has much 
~~ for HD 97048 and 

(80 % of the [OI] 



wider wings (FWHM = 24 and 25 km s~ 
HD 100546), corresponding to a smaller Ri, 
emis sion emanates from the region between 0.8 and 20 AU, Fig. 
7 in lAcke & van denAncker(20Ql))- 



rect im aging of the circumstellar disk (Fig. SI of Lagage et al 
d2006l) ), compared to 160 ±19 °, found bv lAcke & van den Ancker 



(2006) using astrometry. Using the latter value makes the simu- 
lated line profile slightly more double peaked, and thus increases 
the minimal but does not change the best fit parameters. 



6. Discussion 

Fundamental CO emission from Herbig stars, and gaps therein, 
have been detected in more circumstellar disks, e.g. in SR 
21 with an R m of 7.6 ± 0.4 AU for the 12 CO v(l-0) fine 
dPontoppidan et all 120081) . while there is an inner dust com- 
ponent at < 0.5 AU, followed by a dust gap up to 18 AU 
ilBrown et all 120071) . iGoto etalJ (120061) find an inner clearing 
in the CO emission at 11 + 2 AU around HD 141569, and 
suggest the cavity most likely is cleared by photo-evaporation 
working together with viscous accretion. Britt ain et al.l (120071) 
observe the same star, find CO emission starting between 6 and 
15 AU, and explain the emission as due to UV fluorescence. 
Blake & Boogert (2004) study a sample of 5 Herbig stars, not 
known to have disk gaps or inner holes, and find that the CO gas 
traces both collisional driven emission in the inner (< 0.5 - 1 AU) 
disk, as well as CO gas excited by resonance fluorescence (R < 
50-100 AU). These observations suggest that CO depletion in 
the inner disk is associated with large disk gaps or inner holes. 
Finally, in all cases discussed above (except SR 21 for which 
no data exists), the [OI] emission, when detected, arises at sim- 
ilar or higher velocities than the CO fundamental ro-vibrational 
emission. 

Contrary to the above, the lack of CO emission at small radii 
in HD 97048 and HD 100546 can not be caused by a gas-free 
inner disk. [ OI] 63 00 A emission, a tracer of OH molecular gas 



(lAcke et al. , 12005b in the surface layers of flared disks, and in- 
duced by the stellar UV field, has been observed for both stars by 
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Fig. 5. Top The average line profiles of the HD 97048 v = 
1-0, v = 2-1 and v = 3 - 2 transitions, created with 5 
(5 < J < 7), 12 (5 < J < 14) and 7 (14 < J < 21) lines respec- 
tively, are shown with the black histogram. The blue wing of the 
v = 1 — transition is contaminated by telluric CO lines, shown 
with dots. Over plotted is the 6300 A [OI] emission (red dotted 
line), and the model fit to the CO emission with the solid red 
line. Middle The corresponding, continuum corrected (section[3j 
the uncorrected SPP is shown with the dotted line), average as- 
trometric signal. The model is over plotted in red. Bottom The 
average spatial FWHM in black, with the model over plotted in 
red. 




R in [AU] 



Fig. 6. xL- surfaces for the best fitting disk model for, from left 
to right, the line profile, SPP and FWHM of the v = 2 - 1 tran- 
sition of HD 97048. The minimal )c kA for each fit is given in the 
panels, and the 1, 2 and 3 <x lines represent >68%, >90% and 
>99% confidence levels. The best fit parameters for each panel 
are shown by a red circle. To better show the xLa surface land- 
scape we have overplotted the^ ed = 5, 10 (left two panels, ^ d 



1.5, 1.7 in the right panel) contours. 



lAcke & van den Anckerl (12006b . with 80% of the emission orig- 
inating between 0.8 < R < 20 AU. The SEDs of our targets also 
display large NIR excesses, typical for hot dust c l ose to the star 
d Acke & van den Anckerl l2004t iBouwman eTail 120031) . In the 
cases of HD 97048 and HD 100546, the absence of CO emission 
at Small radii can thus not be attributed to the paucity of matter in 
the inner disk, nor to the absence of UV flux. Disk holes and/or 
large gaps therefore are not the only cause of CO depletion in the 
inner disk. Here, we suggest that CO is efficiently destroyed in 
the inner disk. A possible mechanism is the photo dissociation 
of CO and the subsequent format ion of simple organ ic species 
such as C 2 H 2 , HCN and CH 4 (lAgundez etall l2008h . Further 
analysis (e.g. observations of hot transitions of carbon-bearing 
molecules) is needed to determine the reason why the CO gas is 
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Fig. 7. Same as Fig. 0for HD 100546, but the v = 1 - 0, v = 2 - 
1, v = 3 — 2 and v = 4 - 3 profiles are made from 1 1 (0 < J < 7), 
12 (5 < J < 14), 9 (11 < J < 26) and 12 (19 < J < 34) lines 
respectively. 

not detected in the inner regions, and may offer us a glimpse of 
the molecular chemistry of these inner disks. 

The similar line widths and shapes for all vibrational states 
in HD 97048, as well as the similar continuum corrected SPP, 
suggest a common excitation mechanism and location. The best 
fit parameter for the radial dependency of the CO emission in- 
tensity a — 3.0 can guide models that aim to explain the highly 
excited vibrational line emission, such as UV fluorescence. The 
similar line shape, continuum corrected SPP and FWHM of the 
higher (v = 3-2 and v = 4 — 3) vibrational transitions in 
HD 100546 also suggest a common excitation mechanism and 
location, but the v = 2 — 1 and especially v = 1-0 vibra- 
tional transitions deviate significantly. The asymmetry in the 
data on HD 100546, and the poor fit to the model, point to a 
more complicat e d disk geo metry, in agreement with findi ng of 
IBouwman etai] (120031) and I Acke & van den Anckerl d2006l) 
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